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Abstract 



We consider non-perturbative effects at short distances in theories with confinement. 
The analysis is straightforward within the Abelian models in which the confinement arises 
on classical level. In all cases considered (compact U (1) in 3D and 4D, dual Abelian Higgs 
model) there are non-perturbative contributions associated with short distances which are 
due to topological defects. In QCD case, both classical and quantum effects determine 
the role of the topological defects and the theoretical analysis has not been completed so 
far. Generically, the topological defects would result in 1/Q 2 corrections going beyond the 
standard Operator Product Expansion. We review existing data on the power corrections 
and find that the data favor existence of the novel corrections, at least at the mass scale 
of (1-2) GeV. We indicate crucial experiments which could further clarify the situation 
on the phenomenological side. 

Ph! 

1. The perturbative QCD describes basic features of hard processes, i.e. processes charac- 
terized by a large mass scale Q. On the other hand, the perturbative QCD does not encode 
the effects of confinement. Hence there is a growing interest in power corrections to the par- 
ton model which might be sensitive to the nature of confinement (for a recent review see PJ). 
Moreover, in the the Abelian Higgs model the leading power correction at short distances does 
reflect confinement 0. Namely, there exist short strings which are responsible for a stringy 
potential between confined charges at vanishing distances r — > 0. 

In this note we are exploring Abelian models with confinement in a more regular way 
by including into consideration the compact three and four dimensional U(l) theory. In all 
the cases the confinement mechanism can be understood classically. Moreover, we find that 
the confinement results in additional terms in the static potential at short distances. How 
topological defects can be manifested at short distances, is easy to understand on the example 
of the Dirac string. Naively, its energy diverges quadratically in the ultraviolet but in compact 
£7(1) it is normalized to zero changing the power corrections at short distances. As for 
non-Abelian theories, the Dirac strings are also allowed because of the compactness of the 
corresponding U(l) subgroup. In this sense, there is a similarity between non-Abelian and 
Abelian cases. However, there is an important difference as well. On the classical level, the 
Dirac strings may end up with monopoles which have a vanishing non-Abelian action. Thus, 
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the monopoles observed within the U(l) projection of QCD (for a review see, e.g., ||) is a 
result of an interplay between classical and quantum effects. In all the generality, one may say 
that the topological defects in QCD are marked rather by singular potentials than by a large 
non-Abelian action. Singular gauge potentials might be artifact of the gauge fixing and it is 
not a priori clear whether they can result in physical effects. Therefore, we will turn at this 
point to analysis of existing data on the power corrections. The data seem clearly favor the 
novel 1/Q 2 corrections. 

2. First, we will outline very briefly the standard approach to the power corrections which 
allows to account for soft non-perturbative field configurations (for further references see, e.g., 
[|Ij]). Consider first a QED example ||. Namely, let an e + e~ pair be placed at distance r near 
the center of a conducting cage of size L, L 3> r. Then the potential energy of the pair can be 
approximated as 

V e e{r) ps + const ■ L > r (1) 

r L 6 

and the second term is a power correction to the Coulomb interaction. The derivation of ([I]) 
is of course straightforward classically, since the correction is nothing else but interaction of 
the dipole with its images. In the QCD case, one concludes by analogy that the heavy quark 
potential at short distances looks as (for explanations and further references see, e.g. 0): 

JES^QQW = ~ C -y + const -K% CD r\ (2) 

where c_i is calculable perturbatively as a series in a s . Note the absence of a linear correction 
to the potential at short distances. 

On the other hand, Eq. ([!]) can be derived also in terms of one-photon exchange. The power 
correction is related then to a change in modes of the electromagnetic field confined in the cage 
as compared to the case of the infinite space. The change is of order unit at frequencies 
uj ~ 1/L. Similarly, the logic behind Eq. (Q) is that the perturbative gluon propagator is 
modified strongly by at uj ~ ^qcd- m case °f other processes, the relevant Feynman graphs 
can be more complicated of course. The power corrections still correspond to the infrared 
sensitive part of Feynman propagators which are obviously modified by the physics of large 
distances. The Operator Product Expansion (OPE) allows for a regular way to parameterize 
such corrections (for a review see 

3. Intuitively, the power corrections in QCD could be very different from the conducting 
cage case discussed above. Indeed color particles produced at short distances find no cage 
but rather build up the confining field configuration in the course of the interaction between 
themselves and with the vacuum. The complicated space-time picture of interaction in confining 
theories was studied by Gribov JTj. Thus, it could be instructive to analyze the effects of the 
confinement at short distances in some simple models. 

The first example of a theory where the OPE does not work in fact goes back to the paper in 
Ref. @. However, since it has not been discussed in connection with the OPE, we will explain 
this example in some detail. The action is that of free photons: 
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where is the field strength tensor, F^ v = d ll A u — d u A^ 1 . Although the theory looks absolutely 
trivial, it is not the case if one admits Dirac strings into the theory. Naively, the energy 
associated with the Dirac strings is infinite: 

rn 1 f rs TT 2 7 a /magnetic flux \ 2 

^Dirac string = rfrH ~ l ' A I J J ^ 00 ( 4 ) 

where I, A are the length and area of the string, respectively. Since the magnetic flux carried 
by the string is quantized and finite the energy diverges quadratically in the ultraviolet, i.e. in 
the limit A — > 0. However within the lattice regularization the action of the string is in fact 
zero because of the compactness of the U(l), Ref. ||. 

Now, the Dirac strings may end up with monopoles. The action associated with the 
monopoles diverges in ultraviolet, 

£h> ~ i (5) 



8n e 2 a 

where a is a (small) spatial cut off. If the length of a closed monopole trajectory is L, then the 
suppression of such a configuration due to a non-vanishing action is of order 



e ~ exp (—const L/e ). (6) 

On the other hand, there are different ways to organize a loop of length L. This is the entropy 
factor. It is known to grow exponentially with L as ~ exp( const' L ). At some ^crit ^ 1 there 
is a phase transition to the monopole condensation. 



The potential between external test charges is Coulombic at all the distances for e 2 < e 2 



crit 



and linear for e > e crit - Since there are no perturbative graphs at all in the theory with the 
action (|3]) this phenomenon clearly goes beyond the OPE. In this case, however, the violation 
of the OPE is too strong. Indeed, the lattice spacing a is the only dimensional parameter of 
the problem. As a result Coulomb potential is not simply modified by linear corrections but 
rather eliminated for e 2 > e 2 rii at all the distances. 

4. Consider next the Dual Abelian Higgs Model with the action 

5 = /^|^ + i|(9-^)$| 2 + iA(|$| 2 - 57 2 ) 2 }, (7) 

here g is the magnetic charge, F^ u = d fJi A u — d u A fl . The gauge boson and the Higgs are massive, 
my = g 2 n 2 , m 2 H = 2Xn 2 . There is a well known Abrikosov-Nielsen-Olesen (ANO) solution to 
the corresponding equations of motion. The dual ANO string may end up with electric charges. 
As a result, the potential for a test charge-ant i- charge pair grows linearly at large distances: 

V(r) = (W , r — > oo. (8) 

Note that there is a Dirac string resting along the axis of the ANO string connecting monopoles 
and its energy is still normalized to zero. 



3 



An amusing effect occurs if one goes to distances much smaller than the characteristic mass 
scales Ttiy^jj. Then the ANO string is peeled off and one deals with a naked (dual) Dirac string. 
The manifestation of the string is that the Higgs field has to vanish along a line connecting the 
external charges. Otherwise, the energy of the Dirac string would jump to infinity anew. 

As a result of the boundary condition that $ vanishes on a line connecting the charges the 
potential contains a stringy piece at short distances [0: 

limV(r) = — — + cr • r . (9) 

r^O r 

The string tension o"o smoothly depends on the ratio m#/my. In the Bogomol'ny limit (m# = 
my) which is favored by the fits of the lattice simulations [|J the string tension 

(7 « 0"oo, (10) 

i.e. the effective string tension numerically is the same at all distances. 

5. Consider now 3D compact electrodynamics. As is well known ||, the charge-anti-charge 
potential is then linear at large r. Below we consider the string tension cr at small distances 
and show that it has a non-analytical piece. 

As usual, it is convenient to perform the duality transformation, and work with the cor- 
responding Sine-Gordon theory. The expectation value of the Wilson loop in dual variables 
is: 

W = \ I V X e ~ S(X ' VC) > ( U ) 

where 

S ^c) = {£) J d 3 xfyd X ) 2 + m 2 D (l-cos[x-ric))}, (12) 

mo is the Debye mass and S(x, 0) is the action of the model. If static charge and anti-charge 
are placed at the points (—R/2,0) and (.R/2,0) in the xi,X2 plane (X3 is the time axis), then 

Vc = *rctg[^-^] - arctg[^^], -n < Vc < vr. (13) 

Below we present the results of the numerical calculations of the dimensionless string tension, 

a = dE/d(m D R) , (14) 
E = J d 2 x [\0X) 2 + m 2 D (l - cos[ X - %])} • (15) 



Note that the energy E is measured in the units of the dimensional factor (^) 2 (cf. (p~2|) ) . 
Variation of functional fll5D leads to the equation of motion Ax = rn 2 D sin[x — rj c ]. For finite R 
we can solve this equation numerically. The energy E versus m^R is shown on Fig. 1(a). At 
large separations between the charges (rri£>R ^> 1) it tends to the asymptotic linear behavior 
E = 8m£)i? which can be obtained also analytically 0. 

At small distances there is a contributions of Yukawa-type to the energy (|I5D , which should 
be extracted explicitly. Note that in course of rewriting original 3-D compact electrodynamics in 
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Figure 1: (a) The dimensionless string tension (|TJ|) of the charge - anti-charge separated by 
the distance R in 3D compact U(l) theory; (b) The string tension a strm9 (|17]) as a function of 
tudR with corresponding best fitting function (see text). 



the form ([llrlH) the Coulomb potential was already subtracted, so that (|15|) contains Yukawa- 
like piece without singularity at R = 0. It is not difficult to find the corresponding coefficient: 

E = E stnn9 - 2n{K [m D R] + ln[m D R}) (16) 

where K (x) is the modified Bessel function and E strtn9 is the energy of the charge-anti-charge 
pair which is only due to the string formation. The corresponding string tension 

firing = (J + 2 v(-K 1 \m D R} + — (17) 

rri£)K 

is shown on Fig. 1(b). We found that the best fit of numerical data for small values of vtldR is 
by the function a strm9 = const ■ (rri£)R) u which gives v w 0.6. 

Thus the non-analytical potential associated with small distances is softer than in the case 
of the Abelian Higgs model. The source of the non-analyticity is the behavior of the function 
r\c{x\,x-i) eq. (0) which is singular along the line connecting the charges, see Fig. 2(a). 

6. The compact electrodynamics is usually considered as the limit of Georgi-Glashow 
model, when the radius of the 't Hooft - Polyakov monopole tends to zero. For a non-vanishing 
monopole size the problem of evaluating the potential at small distances becomes rather in- 
volved. To avoid unnecessary complications we consider the 3D Georgi-Glashow model in the 
BPS limit. The 't Hooft - Polyakov monopole corresponds then to the fields: 

rr a f 1 



r ytanh(/xr) 

sinh(/xr) J 



A- = -e aic — (-- . £ . ] , A a = 0. (19) 
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(a) (b) 
Figure 2: Function rjc in compact U(l) theory (a) and for extended monopoles (b), eq. ( 22 ) 

The contribution of this monopole to the full non-Abelian Wilson loop W can be calculated 
analytically. If the static charges are placed at points ±R/2 in the {x\,X2) plane the result is: 

W(bi,b2, u) = cos h(fibi) cos h(fib2) + ^-7— y~— sin hiubx) sin/i(/ife 2 ), (20) 

b\ b 2 

here 61,2 = xo ± R/2, bk = |6fe|, xo is the center of the 't Hooft - Polyakov monopole and 

h(x) = / -== — . (21) 

^ ^ -00 V 3:2 + C 2 smn yx 2 + C 2 

One way to represent ( PD[ ) in terms of the function introduced earlier is: 

T)c{x , R, fi) = sign(y) arccos W(b u b 2 , //) . (22) 



In the limit fj,R — > 00 W(foi,&2>A*) cos 77c and r}c(xo, R, fx) coincides with the definition (|13|). 
For small /ii? the function eq. (|22] ) is singular not only between external charges, but also 
outside this region (see Fig.2(b)) although the strength of singularity gets smaller. In the limit 
of vanishing r]c the string tension at small distances <7q apparently goes to zero. 

To summarize, it is natural to expect that in the Georgi-Glashow model the non-analytical 
piece in the potential disappears at distances much smaller than the monopole size. Note, 
however, that to evaluate the potential consistently in this case one should have taken into 
account also the modification of the interaction due to the finite size of the monopoles. 

7. Knowing the physics of the Abelian models outlined above it is easy to argue that the 
perturbative vacuum of QCD is not stable as well. Indeed, let us make the lattice coarser a la 
Wilson until the effective coupling of QCD would reach the value where the phase transition in 
the compact U(l) occurs. Then the QCD perturbative vacuum is unstable against the monopole 
condensation. The actual vacuum state can of course be different from the U(l) case but it 
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cannot remain perturbative. Similar remark with respect to formation of Z 2 vortices was in 
fact made long time ago [IJJ . 



The existence of the infinitely thin topological defects in QCD makes it close akin of the 
Abelian models considered above. However, the non-Abelian nature of the interaction brings 
in an important difference as well. Namely, the topological defects in QCD are marked rather 
by singular potentials than by a large non-Abelian action. Consider first the Dirac string. 

Introduce to this end a potential which is a pure gauge: 

= (23) 

and choose the matrix Q in the form: 

/ cos I sin I e~ iLp \ 

^ = ■ e i9 e ) ( 24 ) 

V — sin I e %H> cos ~ J 

where ip and 9 are azimuthal and polar angles, respectively. Then it is straightforward to check 
that we generated a Dirac string directed along the X3-axis ending at £3 = and carrying the 
color index a = 3. It is quite obvious that such Abelian-like strings are allowed by the lattice 
regularization of the theory. 

The crucial point, however, is that the non-Abelian action associated with the potential 
(p3|) is identical zero. On the other hand, in its Abelian components the potential looks as a 
Dirac monopole, which are known to play important role in the Abelian projection of QCD (for 
a review see, e.g., 0]). Thus, there is a kind of mismatch between short- and large- distance 
pictures. Namely, if one considers the lattice size a — > 0, then the corresponding coupling 
g(a) — ► and the solution with a zero action ( p3|) is strongly favored at short distances. At 
larger distances we are aware of the dominance of the Abelian monopoles which have a non-zero 
action. The end-points of a Dirac string still mark centers of the Abelian monopole. Thus, 
monopoles can be defined as point-like objects topologically in terms of singular potentials, not 
action. 

Similar remarks hold in case of the Z 2 vortices. Namely the Z 2 vortices which have a typical 
size of order Aqq D can be defined topologically in terms of the so called P-vortices which are 
infinitely thin but gauge dependent, see |TI] and references therein. To detect the P-vortices 
one uses the gauge maximizing the sum 

El Tr ^| 2 (25) 
l 

where I runs over all the links on the lattice. The center projection is obtained by replacing 

Ui -> sign (Tr £/,). (26) 

Each plaquette is marked either as (+1) or (—1) depending on the product of the signs assigned 
to the corresponding links. The P- vortex then pierces a plaquette with (-1). Moreover, the 
fraction p of the total number of plaquettes pierced by the P-vortices and of the total number 
of all the plaquettes Nt, was found to obey numerically the scaling law 

V = ^ ~ f(P) (27) 
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where the function f{(3) is such that p scales like the string tension. Assuming independence 
of the piercing for each plaquette one has then for the center-projected Wilson loop W cp : 

W cp =[(l-p)(+l)+p(-l)] A « e~ 2 ^ (28) 

where A is the number of plaquettes in the area stretched on the Wilson loop. Numerically, 
Eq. ( p8|) reproduces the full string tension. 

It is quite obvious that the P-vortices, since they are constructed on negative links, corre- 
spond in the continuum limit to singular gauge potentials of order a -1 . Moreover, the large 
potentials should mostly cancel if the corresponding field-strength tensors are calculated be- 
cause of the asymptotic freedom. The argumentation is essentially the same as outlined above 



for the monopoles, see, e.g. |12| arid references therein. 

At the moment, it is difficult to say a priori whether the topological defects defined in 
terms of singular potentials can be considered as infinitely thin from the physical point of view. 
They might be gauge artifacts. Phenomenologically, using the topologically defined point-like 
monopoles or infinitely thin P-vortices one can measure non-perturbative QQ potential at all 



the distances. It is remarkable therefore that the potentials generated both by monopoles |L3[ 



and P-vortices |TT[ turn to be linear at all the distances measured: 

V non - P ert{r) ~ cr^r at all r (29) 

Note that the Coulomb-like part is totally subtracted out through the use of the topological 
defects. Moreover, the no-change in the slope ( |2"9"D agrees well with the dual Abelian Higgs 
model as discussed above (for alternative approaches see JRl |I5||). 



The numerical observation (^) is by no means trivial. If it were so that only the non- 
Abelian action counts, then the non-perturbative fluctuations labeled by the Dirac strings 
or by P-vortices are bulky (see discussion above) and the corresponding QQ potentials ( |29|) 
should have been quadratic at small r. This happens, for example, in the model with finite 
thickness of Z 2 vortices. Similarly, if the lessons from the Georgi-Glashow model considered 
above apply the finite size of the monopoles would spoil linearity of the potential at short 
distances. 

To summarize, direct measurements of the non-perturbative QQ potential indicate the pres- 
ence of a stringy potential at short distances. The measurements go down to distances of order 
(2 GeV)" 1 . 



8. In view of the results (29) it is interesting to reexamine the power corrections with the 
question in mind, whether there is room for novel stringy corrections. From the dimensional 
considerations alone it is clear that the new corrections are of order <Jq/Q 2 where Q is a large 
generic mass parameter characteristic for problem in hand. Also, the ultraviolet renormalons 
in 4D indicate the same kind of correction, see || and references therein. Note that unlike 
the case of the non-perturbative potential discussed above, other determinations of the power 
corrections ask for a subtraction of the dominating perturbative part and this might make the 
results less definitive. 

(i) The first claim of observation of the non-standard 1/Q 2 corrections was made in Ref. |L7| . 



Namely, it was found that the expectation value of the plaquette minus perturbation theory 
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contribution shows 1/Q 2 behavior. On the other hand, the standard OPE results in a 1/Q 4 
correction. 



(ii) The lattice simulation [18) do not show any change in the slope of the full QQ potential 
as the distances are changed from the largest to the smallest ones where the Coulombic part 
becomes dominant. An explicit subtraction of the perturbative corrections at small distances 



from QQ potential in lattice gluodynamics was performed in ref.[T!|. This procedure gives 
do ~ 50oo at very small distances. 



(in) There exist lattice measurements |2(J of the fine splitting of QQ levels as function of 



the heavy quark mass. The Voloshin-Leutwyler [21] picture predicts a particular pattern of the 



heavy mass dependence of this splitting. Moreover, these predictions are very different from 
the predictions based, say, on the Buchmuller-Tye potential |Z| which adding a linear part to 
the Coulomb potential. The numerical results favor the linear correction to the potential at 
short distances. 

(iv) Analytical studies of the Bethe-Salpeter equation and comparison of the results with 
the charmonium spectrum data favor a non- vanishing linear correction to the potential at short 
distances [031. 



[v) The lattice-measured instanton density as a function of the instanton size p does not 
satisfy the standard OPE predictions that the leading correction is of order p 4 . Instead, the 
leading corrections is in fact quadratic |25 



[w) One of the most interesting manifestations of short strings might be the 1/Q 2 corrections 
to current correlation functions Uj(Q 2 ). It is not possible to calculate the coefficient of front 
of the 1/Q 2 terms from first principles, however, in Ref. [24| it was suggested to simulate 



this correction by a tachyonic gluon mass. On one hand, the tachyonic mass imitates the 
stringy piece in the potential at short distances. On the other hand, it can be used in one- 
loop calculations of the correlation functions. Rather unexpectedly, the use of the tachyonic 
gluon mass {m 2 = —0.5 GeV 2 ) explains well the behavior of Uj(Q 2 ) in various channels. To 
check the model further, it would be very important to perform accurate calculations of various 
correlators Hj(Q 2 ) on the lattice. 

9. As seen from the points (i)-(vi) above, the existence of the novel quadratic corrections is 
strongly supported by the data. There are, however, two caveats to the statement that the novel 
short-distance power corrections have been detected. On the theoretical side, the existence of 
short strings has been proven only within the Abelian Higgs model. As for the QCD itself, the 
analysis is so far inconclusive. On the experimental side, the data always refer to a limited 
range of distances. In particular, the linear non-perturbative potential has been observed at 
distances of order of one lattice spacing which in physical units is about (1 4- 2 GeV) -1 . One 
cannot rule out that at shorter distances the behavior of the non-perturbative power corrections 



changes (see, e.g., [15, |25|). Which would be a remarkable phenomenon by itself. 
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